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ABSTRACT 


Methods  for  the  detection  of  the  optical  signature  of  sub¬ 
marines  by  means  of  laser  techniques  are  discussed,  The 
principal  approaches  considered  are  those  based  on  the  optical 
doppler  effect  and  on  the  optical  matched  filtering.  The 
experimental  results  obtained  in  the  course  of  the  investiga¬ 
tion  are  presented  and  are  discussed  from  the  point  of  view 
of  the  feasibility  of  the  techniques.  It  is  shoim  that,  while 
it  is  difficult  to  extend  the  optical  doppler  method  from 
lab<jratory  to  field  operation,  the  optical  matched-filtering 
method  lends  Itself  to  the  design  of  simple  instrumentation 
for  the  detection  of  the  optical  signature  in  real  time. 
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1*  INTRODUCTION 

The  following  report  covers  work  of  investigation  of  techniques  for 
the  detection  of  submarines  by  means  of  laser  techniques.  The  work  was 
conducted  under  Contract  NOnr  v&b  coiq>leted  in  June  I966. 

The  investigation  has  been  directed  to  the  detection  of  the  "optical 
signature"  of  submarines.  The  equilibrium  configuration  of  the  sea>air 
interface  is  determined  by  the  elastic  waves  propagated  within  the  water 
and  corresponds  to  a  surface  where  the  pressure  equals  that  of  the  atmos¬ 
phere;  thus^  any  perturbation  caused  by  submarines,  vessels,  as  well  as 
winds,  etc.,  contributes  to  the  resultant  configuration.  In  particular, 
a  submarine  in  motion  constitutes  a  distinct  source  of  such  waves  which 
produces  surface  "signatures"  such  as  wakes  and  modifications  of  the  mean 
slopes  of  water  waves. 

Visual  as  well  as  radar  detection  of  submarine  wakes  have  been 
demonstrated;  however,  these  means  of  detection  have  not  proven  suffi¬ 
ciently  practical,  sensitive,  or  reliable.  It  appears  that,  by  recourse 
to  coherent  optical  techniques,  a  greatly  enhanced  sensitivity  may  be 
obtained,  thus  opening  the  way  to  the  design  of  field  equipment  with  suit¬ 
able  operational  characteristics. 

In  the  following,  two  approaches  to  the  problem  of  detection  of  the 
optical  signature  by  meat*  of  lasers  will  be  described*  The  first  is 
based  on  the  measurement  cf  optical  frequency  shifts  produced  by  doppler 
effects,  and  the  second  is  based  on  the  use  of  optical  Intensity  scatter 
as  well  as  optical  matching  filters  and  electronic  systems  of  signal-to- 
noise  enhancement.  He  precede  the  treatment  with  a  review  of  the  basic 
theories  of  hydrodynamics  which  provide  a  quantitative  description  of  the 
Optical-signature  phenomena  at  the  sea  surface. 
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2.  PRINCIPUES  OF  HYDRODYNAMICS 

The  analysis  of  the  sea-surface  waves  has  been  the  subject  of  extensive 
Investigations;  the  basic  theories  have  been  developed  In  hydrodynamics  and 
are  associated  with  the  names  of  Lagrange,  Stokes,  Kelvin,  Lamb,  to  name  only 
a  few.  In  particular,  Kelvin  has  provided  a  beautiful  description  of  the 
physical  phenomena  of  wakes  created  by  ships  in  motion. 

The  equations  of  hydrodynamics  may  be  expressed  either  with  reference 
to  the  motion  of  individual  fluid  particles  (Lagrange  form),  or  with  refer¬ 
ence  to  the  velocity  distribution  at  fixed  points  In  the  medium,  as  functions 
of  time  (Euler  form).  In  general,  Indicating  with  p  ^the  density  of  water, 
with  g  the  pressure,  with  F  the  body  force,  with  a  the  acceleration,  and 
with  V  the  velocity,  one  can  write  the  following  general  equations: 


Equation  of  Motion 

~  grad  p  +  F  -  a 
p 

(1) 

Equation  of  Continuity 

dlv  V  »  0 

(2) 

Equation  of  Irrotational 
Flow 

curl  V  ■  0 

(3) 

As  a  result  of  application  of  Equation  O’Q*  can  express  the  velocity 
vector  as  the  gradient  of  a  single-valued  scalar  potential  function  0  (x, 
y,  *,  t)  ;  l.e.. 


V  <■  grad  0 


(4) 


This  relation  is  very  Important,  because  it  reduces  the  study  of  the  vector 
quantity  v  to  that  of  a  scalar  function;  substituting  Equation  h  into 
Equation  2,  one  has: 

7*0-0  (5) 

i«e.,  the  function  0  is  obtained  as  solution  of  the  linear  scalar  Laplace's 
equation. 

Substitution  of  Equation  h  into  Equation  1  provides  a  representation  of 
the  motion  in  terns  of  0  (Eq.  of  Bernoulli);  for  example,  the  equation  of 
motion  of  a  fluid  particle  subject  to  the  force  of  gravity  is  expressed  as 
follows t 

jf  4  I  (u®  +  V*  +  w*)  +  J  +  gy  -  c(t)  (6) 
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where  u  ,  v  ,  w  are  scalar  velocity  coo^onents,  and  c(t}  is  a  function 
derived  from  time-dependent  external  forces.. 

Thus,  the  solution  of  the  problems  of  hydrodynamics  is  carried  out  deter¬ 
mining  first  the  scalar  potential  function  0  from  Equation  ^  ,  and  secondly 
the  velocity  components  u  ,  v  ,  v  from  Equation  6. 

The  conditions  at  the  boundaries  of  the  sea  medium  are  of  two  types: 
a)  those  corresponding  to  the  sea  bottom,  where  the  surface  is  fixed  and 
therefore  the  normal  consonant  of  the  velocity  is  zero;  these  boundary  con¬ 
ditions  are  expressed  as  follows; 


^  .  0  on  S 


(7) 


b)  those  corresponding  to  the  sea-air  Interface,  where  the  vertical  component 
of  the  acceleration  is  zero,  because  the  surface  is  ''free”  and  assumes  a  form 
such  that  a  ■  0  .  Referring  to  a  Cartesian  reference  system  of  axes,  with 
origin  at  the  undisturbed  water  surface  and  y-axis  directed  vertically  upwards, 
the  latter  boundary  condition  is  expressed  as  follows; 


^x  \  '  ^y  ^t  “  ®  ® 

where  0,  0,  0,  ij*,  n  ,  n  are  partial  derivatives  with  respect 

to  the  variables  x  ,  y  t  ,  *t  ,  and 


y  -  ii(x  ,  *  ,  t)  (9) 


is  the  equation  of  the  "free^  surface  S  . 

The  above  relationships  provide  the  mathematical  means  for  the  determina¬ 
tion  of  the  propagation  of  traves  and  of  the  configuration  of  the  sea  surface. 
For  exaople,  progressive  surface  gravity  waves,  where  the  gravity  force  is 
predominant  ever  the  surface  tension,  are  represented  with  the  following  velo¬ 
city  potential  function; 


0  •  A  cos  h  p(yfh)  cos  (^x  l  ot  •f  o)  (10) 

I 

where  h  is  the  d^th  of  the  sea,  a  is  a  phase  angle,  and  the  quantities 
P  and  0  have  the  following  significance; 
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0  ■  gP  taa  h  ph 


If  the  depth  h  is  very  large. 


The  velocity  component  along  the  x  direction  is  expressed  as  follotfs: 


V  B  u  «>  r- 
P 


a-  Vi 


More  generally,  if  the  surface  tension  T  is  considered,  the  latter  rela¬ 
tion  is  expressed  as  folloi*rs: 


The  function  14  is  represented  graphically  in  Fig.  1  versus  the  wave¬ 
length  X  ;  it  is  noted  that  the  velocity  changes  markedly  depending  on 
whether  the  surface  tension  T  or  the  gravity  g  is  the  restoring  force. 

The  former  case  corresponds  to  capillary  waves,  and  the  latter  to  gravity 
waves. 

The  frequency  of  oscillation  of  the  waves  is  computed  with  the  familiar 
relationship  Xf  ■  v  ,  and  thus  varies  with  the  wavelength,  with  the  follow¬ 
ing  functional  relationship: 


*  "yTxT' 

With  reference  to  Fig.  1|  it  Is  noted  that  a  minimum  velocity  of  20-30  cm/ 
sec  occurs  in  correspondence  of  X  ■  I  r-  2  cm  ;  capillary  waves  display 
wavelength  values  less  than  1  cm  ,  and  frequency  larger  than  10  ~  13  c/sec. 
Schooley  has  pointed  out^'^^3  that,  while  small  waves  (both  capillary  and 
gravity  type)  have  epproximately  sinusoidal  profiles,  large  waves  of  gravity 
type  have  a  modified  trochoidal  shape,  sharper  near  the  crests  and  flatter  in 
the  troughs,  while  large  waves  of  capillary  type  have  sharp  troughs  and  wide 
crests  and,  in  general,  larger  average  slope  than  gravity  waves  (Fig.  2).  The 
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surface  slope  also  depends  upon  the  wind  and  statistically  is  peaked  and 
skewed  in  the  down^wind  direction  (compared  to  a  Gaussian  distribution  - 
Fig.  3). 

As  far  as  the  wake  pattern  is  concerned,  it  is  sufficient  to  summarize 
the  results  derived  from  Kelvin^  .  Indicating  with  C  the  course  of  a 
ship,  and  assming  that  the  origin  of  the  coordinates  axes  moves  with  the 
ship  (Fig.  h),  the  curve  C  may  be  represented  parametrically  with  the 
equations : 


V 


Xj  «  x^(f) 

O^t^T 

*1  “  *1^-*^^ 


(16) 


where  t  is  the  time  variable.  Kelvin  has  shown  that  the  vertical  profile 
of  the  sea  surface  y  “  tj  (z,  t)  may  be  expressed  as  an  integral; 

3  2 

n  (*i  t)  -  k^J  ^  sin  dt  (17) 

/o  ' 

where  r  yields  the  distance  between  a  point  (x,z)  and  a  point  (x^^,  z^^) 
on  the  course; l.e.. 


The  Integral  (17)  was  computed  approximately  by  the  method  of  "station¬ 
ary  phase".  Loci  of  constant  ph;ise,  for  a  number  of  phase  values  differing 
by  2  IT,  are  slwwn  In  Fig.  5.  for  the  case  of  a  straight  course  at  constant 
spaed  and,  in  Fig.  6,  for  the  case  of  a  circular  course  at  constant  speed.  In 
the  former  example  (Pig.  5)«  the  curves  are  represented  in  terms  of  the  angle 
9  ,(^ig*  4),  with  the  following  equation: 


X  -  ^  (2  cos  6  -  cos^  6) 
a  2  ' 

a  «  I  cos  6  sin  6 


(19) 


where  a  is  a  parameter  depending  upon  the  velocity  «sd  having  the  dltaeoalons 
of  length.  These  curves  constitute  two  systems  of  waves,  respectively  diver¬ 
gent  and  transverse,  which  cover  a  V-shi^ed  region  celled  the  "influence  re¬ 
gion*,  and  join  at  the  boundary  of  such  region;  ttM  points  of  Junction  are 
cusps  and  correspond  to  the  condition  9  ■  9*  ,  where 


( 
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3  8ln^  e*  -  I  *  (20) 

It  can  be  shown  that  In  first  approximation  the  disturbance  at  a  point 
F(x,y)  is  determined  only  by  certain  discrete  points  Q.  on  the  ship  course 
C  given  by  the  relation: 

r  •  ^  c  t  cos  6  (21) 


where  c  and  t  are  respectively  speed  and  time;  In  fact,  the  contributions 
from  all  other  points  cancel  out  because  of  mutual  interference.  The  points 
Q  are  called  "influence  points";  for  a  given  point  Q  ,  one  can  compute  from 
(21)  the  values  of  the  coordinates  (r,  8)  of  points  P  for  which  Q  is  the 
influence  point.  Such  points  £  are  found  on  a  circle  (Fig.  6)  with  diameter 
tangent  to  the  course  C  at  Q  ,  where  Q  is  on  the  circumference.  This 
property  permits  one  to  construct  geometrically  the  outline  of  the  region  in 
which  the  disturbance  caused  by  the  ship's  motion  is  confined.  For  example, 
in  Fig.  7i  the  ease  of  a  straight^path  course  is  shewn;  circles  of  radius  1/k 
c  t  are  drawn  at  various  points  Q  and  delimit  by  envelope  the  "wake"  bound¬ 
aries.  He  see  that,  if  the  speed  c  is  constant,  th*  boundaries  are  two 
straight  lines  with  semiangle  t  19*28'  independent  of  the  speed  c  .  The 
importance  of  this  result  for  application  to  detection  is  obvious. 


It  can  be  shown  that  the  curves  of  constant  phase  of  a  wake  are  ortho¬ 
gonal  to  the  lines  drawn  back  to  the  influence  points.  This  is  shown  graphi¬ 
cally  in  Pig.  8  and  is  demonstrated  by  differentiation  of  Equation  19*  One 
finds: 


and 


^  “  -  I  (3  8-1)  sin  8 

^  "  5  (3  sin®  8-1)  cos  8 

de 


da  _ _ l_ 

dx  "  ’  tan  8 


(22) 


(S3) 


From  Equations  22  and  23,  it  is  also  found  that  the  points  8*  defined  by 
Equation  20  are  singular  and  correspond  to  the  cases  in  which  the  influence 
points  of  the  divergent  and  transversal  constant-phase  waves  coincide. 

Finally,  the  amplitudes  of  the  wake  waves  can  be  confuted  approximately, 
although  laboriously.  In  ^neral,  it  is  found  that  they  are  Inversely  proper- 
tiMial  to  the  square  root  of  the  distance  from  the  ship. 
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The  doppler  effect  may  be  used  to  measure  the  relative  velocity  of  a 
target  with  respect  to  observer  as  a  function  of  the  shift  of  the  frequency 
of  a  reflected  electromagnetic  wave  as  seen  by  an  observer.  Assuming  that 
the  component  of  velocity  along  the  path  of  propagation  and  reflection  is  v, 
the  frequency  f  »of  electromagnetic  wave  scattered  from  the  target  is: 


Expanding  in  series  and  neglecting  higher-order  terms,  one  has: 

f  -  f  (1  -  J) 

or  (25) 

£*  -  f  ■  f^  «  £  v/c 


More  generally,  assuming  that  the  angle  between  the  velocity  vector  v  and 
the  propagation  path  Is  o  ,  one  has: 


where  X  is  the  wavelength  'o£  the  wave. 

In  practice,  the  scatter  phenomena  are  accompanied  by  more-complex  phase 
shi£t8  depending  upon  the  path-length  spread  due  to  beam  widening  and  to  time- 
varying  o£  the  distance.  Assuming  that  the  width  of  the  received  beam  is  @ 
(Fig.  9)f  the  distance  R  varies  between  two  limits  R^  and  Ry^  so  that,  the 
corresponding  phase  of  the  returned  wave  is  also  varied  between 


(87) 


la  additlm,  th.  dlstue.  R  VMl..  ia  .  WMll  tin.  iattrval  with  th. 
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Thus,  If  the  transmitted  wave  is  expressed  as 
m 


“T 


-  sin  2irf 


(29) 


n»0 


the  received  wave  assumes  the  expression 


®  ItTT 

G  sin  2Trf  (t-t^)  -  ^  (t^)  -  v(t-t^)  cos 


(30) 


where  6  is  a  spread  loss  factor  depending  upon  target-scattering  properties. 
The  latter  equation  shows  that  the  reflected  beam  is  characterized  by  a  fre¬ 
quency  spectrum  rather  than  by  a  single  frequency;  in  fact,  writing  Eq.  30 
as  follows: 


one  observes  that  a  spread  in  frequency  and  in  phase  is  obtained,  depending 
V  '  upon  the  magnitude  of  the  beamwidth  6  .  Reducing  the  latter  by  recourse  to 
a  laser  beam  is  obviously  beneficial. 

There  are  two  methods  which  may  be  used  to  apply  laser  doppler  shift 
techniques;  one  consists  in  utilizing  an  unmodulated  laser  carrier  and  in  meas¬ 
uring  the  shift  of  the  optical  frequency  due  to  doppler  effect,  and  the  other 
one  consists  in  utilizing  a  laser  beam  modulated  with  a  microwave  frequency  and 
In  measuring  the  doppler  shift  of  the  modulation  rather  than  that  of  the  car¬ 
rier.  The  first  technique  requires  recourse  to  optical  heterodyning  and  takes 
advantage  of  the  \ery  high  value  of  the  optical  frequency;  for  example,  in  the 
case  of  a  HeNe  laser,  one  has: 


m 


G  E  sin 
n 


4 

i 


2V  cos 


2ir  (f  + 


(‘■‘o)  - 


kv 


t'o) 


(31) 


X  -  6328  A  , 


f  -  4.75  X  10^^  c/s 


The  second  technique  requires  recourse  to  microwave  modulation  and  detection 
of  laser  carriers,  but  may  utilize  non-heterodyne-type  detection;  the  basic 
frequency  Involved  In  the  doppler  effect  Is  that  of  the  modulating  signal, 
which  is,  of  course,  many  orders  of  magnitude  lower  than  that  of  the  laser 
carrier.  Thtis,  the  laser  beam  in  this  case  may  be  useful  as  a  narrow-beam 
carrier,  since  It  permits  the  utilization  of  high-gain  transmitting  and  receiv¬ 
ing  optical  antennae.  There  appears  to  be  considerable  merit  in  the  use  of  the 
second  approach,  mainly  because  it  does  not  require  heterodyne  detection;  how¬ 
ever,  difficulties  are  found  in  the  design  of  efficient  microwave  modulators 
and  detectors. 
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The  first  approach  (i.e.,  unmodulated  laser  carrier  and  direct  optical 
heterodyning)  was  investigated  in  our  experimentation. 

A  block  diagram  of  an  optical  laser  radar  is  shown  in  Fig.  10;  it  con¬ 
sists  of  a  CW  laser,  of  an  optical  frequency  translator,  of  an  optical  sum¬ 
ming  device  such  as  an  interferometer  in  which  the  reflected  radiation  is 
added  in  proper  phase  with  the  "local  oscillator"  component,  and  of  a  photo¬ 
detector  and  radio-frequency  receiver. 

Assuming  that  the  laser  beam  is  very  narrow,  the  signal  returned  from 
the  wave  surface  of  the  sea  will  be  deviated  periodically  following  the  con¬ 
tinuous  change  of  the  orientation  of  the  surface  of  incidence  (Fig.  11) ;  if 
the  wave  surface  moves  from  Position  A  to  Position  B,  the  reflected  laser 
beam  is  deviated  from  direction  r^  to  direction  r^  .  Therefore,  the 
returned  light  enters  the  receiver  aperture  at  different  points  and  may  even 
illuminate’  it  only  periodically  for  finite  time  intervals.  This  angular 
deviation  complicates  the  heterodyne-detection  phenomenon  and  must  be  taken 
into  account;  in  the  case  of  wlde-beam  return,  the  receiver  aperture  is  con¬ 
tinuously  illuminated,  but  the  angle  of  arrival  may  vary  between  values  too 
large  to  satisfy  the  conditions  for  heterodyne  detection,  so  that  the  detected 
signal  is  a  pulse  whose  spectrum  depends  upon  the  amplitude  and  period  of  the 
sea-surface  waves.  A  pulse  is  obtained  when  the  laser  beam  has  very  small 
width,  and  the  receiver  aperture  is  Illuminated  periodically. 

A  rough  analysis  of  this  phenomenon  can  be  developed  assuming  that  the 
surface  waves  have  a  triangular  form  and  that  reflected  laser  beam  is  very 
narrow  (Fig.  12);  indicating  with  h  the  peak-to-peak  amplitude  and  with  A 
the  wavelength,  the  slope  of  the  wave  is  given  by  the  relation: 

tana  «  ~  (32) 


In  order  that  a  return  is  obtained  in  the  direction  of  the  incident  beam,  the 
latter  must  be  directed  within  the  angle  characterized  by  the  vertical  and 
the  normal-to-the-wave  slope;  the  maximum  angular  deviation  between  reflected 
and  incident  beam  is  2a  ,  and  the  corresponding  angular  velocity  of  the 
ref lected  beam  Is; 

o  -  to/T 

if  the  sea  waves  are  stationary  and  have  period  T  •  The  duration  of  the 
detected  pulse  Is  determined  by  the  ratio  between  the  angular  aperture  of  the 
photodetector  entrance  pupil  and  the  angular  velocity  of  the  reflected  beam; 

l.a.. 


t  -  D/R  -  DT/llOR  (33) 

(' 


-9- 


CONHDENTIAL 


CONFIDENTIAL 


where  D  is  the  dlsmeter  of  the  entrance  pupil  and  R  is  the  distance  of 
the  same  from  the  point  of  incidence  on  the  sea  surface.  In  practice^  the 
duration  of  the  detected  pulse  depends  also  upon  the  width  of  the  laser  beam 
and  upon  the  velocity  of  propagation  of  sea  waves. 

Recapitulating,  although  the  transmitted  beam  is  continuous,  the  reflected 
beam  gives  rise  to  a  pulsed  detected  signal,  whose  characteristics  of  periodi¬ 
city  and  duration  depend  upon  the  oscillatory  motion  and  translational  velocity 
of  the  sea  waves,  as  well  as  on  the  width  and  relative  orientation  of  incident 
and  of  scattere*  beam.  Further  modifications  must  be  considered  when  the  laser 
beam  is  made  to  scan  a  given  area  of  the  sea  surface. 

Following  the  concepts  of  design  of  the  doppler  laser  shown  in  Fig*  10, 
the  received  pulse  of  duration  T  and  frequency  f^^  +  f^  is  mixed  with  the 
local-oscillator  signal,  which  has  either  frequency  f^  or  f^  +  f^  ;  in  the 
first  case,  one  has  a  homodyne  detector,  and  in  the  second  case,  one  has  a 
heterodyne  detector,  based  on  the  use  of  a  "frequency  translation"  of  the  laser 
carrier.  j 

In  general,  it  is  desirable  that  the  laser  output  be  single  mode  and 
single  frequency;  single-mode  operation  (generally  TEMoqj^)  possesses  the  high 
est  possible  degree  of  coherence,  both  spatial  and  temporal,  and  single- 
frequency  operation  provides  unequivocal  output  of  the  heterodyne  detector. 

Single-mode  operation  may  be  obtained  by  proper  design  of  the  interfero¬ 
metric  cavity^  or  by  application  of  certain  methods  of  frequency- locking  and 
cancellation^^ '  ;  single  frequency.and  frequency  stabilization  means  also 
have  been  investigated  extensively. 

Assuming  that  the  relative  radial  velocity  of  the  reflecting  waves  is  v  , 
the  doppler  frequency  shift  which  appears  at  the  output  of  the  photodetector 
is: 


-  2v/x  ;  (34) 

for  example,  if  v  ■  Im/sec  and  X  ■>  0.63p  ,  one  has; 
fj  -  3«17  Bc/sec 

The  receiver  which  follows  the  photodetector  must  be  able  to  operate  in  a 
frequency  range  encompassing  the  expected  doppler  frequencies;  if  a  frequency 
translation  of  the  local  oscillator  is  used,  the  receiver  frequency  range  is 
reduced  to  f^  «  f  ^  .  Such  frequency  translation  may  be  obtained  by  modulat¬ 
ing  the  reflecting  mirror  of  the  IWy^-Green  Interferometer  used  to  obtain 
superposition  between  reflected  and  local  oscillator  beams  (Fig.  I3);  the 
bandwidth  of  the  receiver  must  be  sufficient .  to  accept  themergy  of  the 
returned  pulse,  which  has  duration  t  ;  i,e.,  Af  ■  1/t.  There  follows  that 
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the  noise  energy  is  determined  by  ^  f  ,  The  receiver  may  be  either  a  super¬ 
heterodyne  with  a  scanning  local-oscillator  frequency,  or  a  broadband  type 
with  channelized  input. 


The  operation  of  the  optical  heterodyne  detector  depends  upon  the  coher¬ 
ent  nature  of  the  incident  beams.  Assuming  that  the  signal  and  the  local 
oscillator  waves  are  cophasal  over  the  entire  surface  of  the  photodetector 
and  that  their  amplitudes  and  frequencies  are  respectively  ,  eOg  ,  Pj^  , 
0)  X  '  finds  that  the  input  current  resulting  from  mixing  and  optical 
detection  is; 


(35) 


kT 


1 


**s  +  ^r  ^  (  “s  " 


The  ac  component  of  this  current  has  mean  square  value: 


2(5^)' 

^  hv  ^ 


P  P 
8  r 


(36) 


The  mean  square  noise  currents  due  to  signal,  local  oscillator,  dark 
current,  radiation  background,  etc.,  are  expressed  as  follows: 

^8^  +  +  ^d^  +  “  2  q  (ig  +  ly  +  +  1^)  ££  (37) 

where  Af  is  the  receiver  bandwidth,  and  I  ,  I  >  Ij  >  Iw  the  detected 

currents  due  respectively  to  signal,  local  osSillaior,  thermionic  emission, 
background  radiation. 


The  quantity  q  is  the  quantum  efficiency  of  the  photpdetector;  for  Si 
and  for  surfaces,  one  has  respectively  q  «  1.5  x  10"^  and  k  x  lo”®  , 
The  signal-to-nolse  ratio  resulting  from  the  consideration  of  the  abo\e  rela¬ 
tion  is; 


_ 

2,  [l.  +  +  Ij  +  H] 


If  the  local  oscillator  power  is  sufficiently  large,  the  above  relations 
BMy  be  sii^plified,  letting: 
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I  •  ^  P 
*r  hv 


(39) 


_£  '  8 
N  "*  hvAf 


This  relation  shows  that  the  signal'to-nolse  ratio  of  coherent  detec¬ 
tors  is  directly  proportional  to  the  quantum  efficiency  of  the  detector;  if 
the  local  oscillator  power  is  sufficiently  large,  the  slgnal-to-noise  ratio 
is  independent  of  all  noise  sources  except  the  photon  noise  hv  Af  .  The 
latter  result  is  important  and  represents  one  of  the  major  advantages  of 
coherent  versus  noncoherent  detectors* 

However,  in  order  to  achieve  operation  of  coherent  detectors,  it  is 
necessary  to  assure  very  strict  requirements  of  optical  alignment.  For 
example,  if  the  target  region  illuminated  by  the  laser  beam  has  linear  dimen¬ 
sion  of  the  order  d  and  is  at  a  distance  R  ,  then  a  necessary  condition 
for  the  realization  of  spatial  coherence  is  that  the  receiver  aperture  have 
diameter  D  where 


D,  <  XR/d 

V 


If  the  transmitter  aperture  diameter  is  D  and  the  beam  is  collimated, 
the  diameter  of  the  spot  size  illuminated  on  the  target  is: 


d  -  +  Xr/D^- 


Bence,  the  diameter  of  the  receiver  aperture  must  satisfy  the  condition: 

D,. 

D  <  t - ^ -  (*>2) 

p 

At  large  ranges  (i*e.,  when  R»  /X  )  #  one  finds  the  simplified  result; 


For  example,  if  the  transmitted  laser  beam  has  diameter  0. 5  cm  ,  the  receiver 
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aperture  diameter  must  be  smaller  than  0. ^  cm  .  A  practical  system  for 
optical  heterodyne  doppler  measurements  is  shown  in  Fig.  1^. 

A  study  of  the  directional  sensitivity  of  the  receiver  may  be  made, 
assuming  that  the  transmitted  and  received  waves  are  plane,  with  propagation 
vectors  and  kg  ;  if  the  detector  surface  is  placed  at  Z  «  0  (Fig.  1?]^ 
the  defected  current  possesses  the  following  expression: 

Dr 

sin  b  ^ 

i(t)  «  cos  (ffij  -  fl:^)t  +  — g -  (44) 

F* 

where 


b  «»  -  kj  sin  0j^  +  kg  sin  Sg  (45) 

i 

A  plot  of  the  amplitude  of  i(t)  versus  b  D^/2  is  shown  in  Fig.  l6; 
it  is  seen  that  the  amplitude  becomes  zero  periodically  ift  correspondence  of 
the  values  b  *•  2n  ir  .  If  one  assumes  that  6^^  »  0  (i.e.,  that  the  local 
oscillator  beao^ is  aligned  so  that  its  direction  of  propagation  is  along  the 
normal-to>the»reeeiver  surface),  the  first  zero  of  the  amplitude  of  i(t) 
occurs  when  kg  sin  Sg  »  kg  ©g  ■  2  ir;  i.e,, 

eg  i  x  /D,  (46) 

,  *4 

For  Instance,  letting  X  O.op  ,  D  ■  0. 5  cm  ,  one  has  Og  -  1.2  x  10 

rad,  A  typical  plot  of  the  angular  sensitivity  of  the  detector  given  in 
Fig.  l6;  this  shows  that  the  heterodyne  detection  of  two  coherent  plane  wave 
beams  Is  extremely  sensitive  to  angular  orientation. 

In  general,  an  heterodyne  receiver  is  equivalent  to  an  antenna;  Indlcst* 
ing  with  A^  (fi)  the  effective  aperture  or  capture  cross-section  for  plane 
waves  Incident  from  direction  0  ,  the  requirement  of  wavefront  alignment  be¬ 
tween  signal  and  local  oscillator  beams  is  translated  into  the  condition 

While  it  is  possible  to  trade  off  effective  aperture  and  angular  direc¬ 
tion  of  arrival,  it  Is  not  possible  to  modify  the  magnitude  of  the  product 
which  is  very  small.  The  condition  (Eq.  47)  imposes  a  severe  limitation  on 
the  practical  utilisation  of  optical  heterodyne  detection  for  doppler  laser 
radars;  this  is  especially  true  when  the  distances  R  involved  are  larger 
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than  ^  •  100  ft  and  when  the  turbulence  of  the  atmospheric  path  and  the 
randomness  of  the  reflecting  surfaces  reduce  the  spatial  coherence  of  the 
laser  beam. 

One  additional  limitation  of  the  design  Is  provided  by  the  time-dependent 
characteristics  of  the  received  pulse  and  by  the  characteristics  of  the 
doppler  signal.  In  fact.  If  the  pulse  duration  is  T  ,  the  minimum  doppler 
frequency  which  can  be  detected  is  1/t  ;  this  represents  the  speed  resolu¬ 
tion  of  the  radar.  For  example.  If  T  *■  1  pse'c  and  X  »  0.63  i  the  minimum 
doppler  frequency  is  1  mc/s  and  the  corpbsponding  speed  resolution  Is  ^ 
0.315  m/sec  B  1  .  Since  the  duration  of  the  reflected  received  pulse 

signal  Is  determined  by  the  oscillatory  and  translational  motion  of  the  sea 
waves,  the  limit  of  resolution  of  the  radar  cannot  be  defined  clearly. 
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k,  EXPERIMENTAL  INVESTIGATION  OF  THE  DOPPLER  FREQUENCY  MEASUREMENTS  BY 

OPTICAL  HETERODYNE  DETECTION 

The  experimental  set'up  built  for  the 'investigation  of  the  doppler  laser 
ASW  T^adar  is  shown  in  Fig.  17*  It  consists  of  a  HeNe  laser  operating  at 
X  =  0,6328^*;  the  optical  cavity  length  is  1.20  m  ,  and  the  mirrors  are 
hemisphcridal  (one  flat  and  the  other  with  radius  2  m  ).  The  diameter  and 
divergence  of  the  collimated  beam  are  respectively  0.5  cm  and  2  x  lO”  rad. 
The  beam  enters  a  Twyman-Green  interferometer  consisting  of  a  beam  splitter 
having.  50'^.  ^ransmisslvity,  and  is  reflected  back  at  the  end  of  two  arms  by 
means  of  d,ielectric  coated  flat  mirrors.  The  latter  are  mounted  on  piezo¬ 
electric  crystals  and  are  positioned  by  means  of  a  3*niicrometers  control  sys¬ 
tem.  A!  photomultiplier  with  Si  surface  has  been  used  as  a  detector;  this 
was  provided  with  suitable  optical  filter  and  iris. 

Experiments  were  conducted  with  the  interferometer  arms  of  equal  length 
(18  cm)  f  and  with  one  arm  of  length  up  to  51"  &nd  the  other  of  length 
18  cm  .  In  addition,  reflections  from  a  water  tank  were  obtained. 

Studies  of  heterodyne  detection  with  equal  arm  length  interferometer  were 
made  In  the  homodyne  mode  and  in  a  heterodyne  mode  obtained  by  means  of  fre¬ 
quency  translation  of  the  "local  oscillator"  wave.  In  the  first  case,  the 
"local  oscillator"  mirror  was  maintained  at  rest,  and  the  "target"  beam  was 
made  to  vibrate  sinusoidally,  using  the  piezo. lectrlc  crystal  control.  The 
frequency  of  oscillation  was  varied  in  the  range  0  -  1  me  .  After  fringes 
of  interference  were  obtained  on  the  surface  of  the  photodetector  in  the 
static  condition,  clear  mixing  was  seen  to  occur.  The  photomultiplier  output 
examined  with  a  spectrum  analyzer  indicated  the  presence  of  Fourier  spectra 
consistent  with  the  sinusoidal  phase  modulation  realized  with  the  vibration 
of  the  mirror. 

From  theory,  one  has  that  the  analytical  representation  of  the  phase- 
modulated  light  wave  is; 


a  ■  A- sin  (a>  t  +  iJ-  sina)  t) 

*  O  01  n 


m 


Expanding  in  a  Fourier  series,  ora  has; 


a  -  ^  ®o  *  L  ^o  ^^m^  *  ^  ^  ®bi  *  '*■  •** 

•f  cos  t  2  +  2  ^m  ^ 

i.a.,  using  known  trigonooetric  expMslons: 


I 


( 
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+  *^3  *\  I  ®^"  K  3a>Jt  -  sin  (to^  -  3coJtl 

•f  . . . 


The  spectrum  resultlnj;  from  this  type  of  modulation  consists  of  sidebands  at 
frequencies  (q  ±  not  >  whose  amplitudes  depend  upon  the  value  of  0  ;  an 

cxan^le  corresponding  to  0^  *>  10  Is  shown  In  Fig.  18. 

Experiments  with  frequency  translation  of  the  local  oscillator  waves  were 
conducted  by  modulating  the  mirror  of  the  "local  oscillator"  arm  with  its 
piezoelectric  mirror.  The  application  of  a  sinusoidal  modulation  again  results 
in  a  spectrum  of  the  type  shewn  above;  in  order  to  avoid  spreading  the  energy 
cn  a  large  number  of  sidebands,  it  is  desirable  to  use  either  a  small  value 
of  0  or  a  suitable  waveform  of  the  modulation  wave;  for  example,  a  ramp¬ 
like  waveform  would  produce  a  single  frequency,  a  sawtooth  waveform  would  pro* 
duce  few  spectral  components,  etc. 

Experimental  verifications  of  the  above  result  are  shown  in  Fig.  19;  20, 
and  21,  where  a  modulating  frequency  of  100  kc  was  used  for  the  "local  oscil¬ 
lator"  wave  and  a  modulating  frequency  of  10  kc  or  9  he  was  used  for  the 
"target"  wave.  The  spectrum  analyzer  was  Panoramic  Spectrum  Analyzer  >!d  SPA 
3125,  with  a  frequency  range  0  -  23  me /s  .  Pigs.  I9,  20,  and  25  correspond 
to  heterodyne  mixing  with  100  kc/'s  modulation  of  the  local  oscillator  and 
5  kc/s  or  10  kc/s  modulation  of  the  "target"  mirror.  Note  that,  in  Fig. 

21,  the  magnitude  of  0  is  selected  so  that  the  carrier  anq>lltude  of  the 
sideband  becomes  zero;  ?n  Fig.  20,  a  single  sideband  at  100  kc  ,  with  the  fre¬ 
quencies  resulting  from  beating  with  it,  is  shown. 

Mixing  of  the  two  incident  waves  was  obtained  also  with  greatly  elongated 
arm  of  the  "target"  wave;  in  particular,  a  one-way  distance  of  5  <1^  vas 
realized.  At  large  distances,  the  alignment  of  the  Interferometer  becomes  more 
critical,  and  the  fringe  separation  decreases  with  constant  optical  magnifi¬ 
cation. 

The  experimental  set-up  for  reflection  from  water  is  seen  in  Fig*  22.  In 
order  to  permit  incidunce  at  sp  angle  for  which  a  reflected  beam  direction 
coincidental  with  that  of  the  transmitted  beaat  may  be  obtained,  the  "depres¬ 
sion"  angle  had  to  be  of  the  order  of  0-10*  •  This  angle  is  shown  in 
Fig*  22,  where  a  practical  set-up  involving  the  use  of  an  intermediate  mirror 
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is  indicated.  Reflections  from  the  surface  of  the  water  were  obtained,  and 
heterodyning  was  demonstrated  easily  in  the  absence  of  water  waves. 

Water  waves  of  stationary  and  of  traveling  type  were  produced  by  means 
of  a  mechanical  stirrer;  frequencies  of  the  order  of  5  **  10  cycles  were 
obtained.  In  addition,  oscillations  of  much  higher  frequency,  though  of  much 
smaller  anq>litude,  were  produced  by  means  of  piezoelectric  transducer  action. 
The  heterodyne  detection  of  the  beams  reflected  from  the  water  waves  was  not 
too  satisfactory.  The  output  of  the  multiplier  appeared  pulsed,  as  discussed 
in  the  previous  section;  the  duration  of  these  pulses  was  a  function  of  the 
geometry  of  the  system,  but  was  too  short  to  permit  detection  of  the  doppler 
signal. 
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5.  .DETECTION  OF  SIGNATURE  BY  OPTICAL  FILTERING  TECHNIQUES 

Because  of  the  limitations  encountered  in  the  application  of  doppler 
frequency  measurements,  a  technique  based  on  the  detection  of  the  optical 
signature  by  optical-filtering  techniques  was  investigated.  During  the 
period  covered  by  this  report,  the  technique  was  investigated  theoretically; 
equipment  was  built  to  provide  a  TV-llke  presentation  of  the  sea-surface 
waves  by  Illumination  and  scatter  with  a  scanning  laser  beam.  Experimental 
studies  of  matched  optical  filtering  and  of  electronic  filtering  should  be 
conducted  in  the  following  period. 

The  concepts  of  optical  spatial  filtering  are  based  on  the  transform 
properties  of  lenses;  ideally,  a  lens  system  transforms  the  light  amplitude 
distribution  in  the  front  and  back  focal  planes  by  means  of  a  Fourier  inte¬ 
gral  relationship.  The  basic  observation  of  this  property  and  its  applica¬ 
tion  to  two-dimensional  filtering  of  signal  patterns  was  first  Introduced  by 
Duffieux  in  France. (1946)  and  was  developed  successively  by  Mfurechal^  , 
0*Neill^®  ,  Gabor ,  and  many  others. 

Consider  the  system  of  Fig.  23>  which  consists  of  an  ideal  lens  of 
aperture  D  ,  illuminated  by  a  monochromatic  unpolarized  wave  having  a  cross- 
sectional  complex  modulation  E^f  (x^y)  ,  at  the  back  focal  plane;  i.e., 

(**  yf  t)  •  Bp  *0  (x>  y)  exp  J  (mt  +  0  (x,y)  )  (51) 


At  the  front  focal  plane,  the  light  distribution  is  expressed  as  follows; 

B,  (*'  y')  -  B^  jj  f^  (x,y)  exp  J  (^  (x,y)  +  x+o^  y)  )  dx  dy 

'  (5S) 

where  ■  Sir  x'A  B  *  K,  ■  «^e  called  "radian  spatial  frequencies 

and  x*  ,  y*  are  coordinates  in  the  front  focal  plane.  The  above  integral 
expression  (52)  is  extended  over  the  pupil  P  of  the  system  and,  except  for 
the  latter  limitation  represents  the  Fourier  transform  of  the  function 
^  (x»y)  •  Thus,  normalizing  the  Integral  (52),  one  has; 


This  eonstitutes  a  two-dimans ional  display  which  may  be  retransformed  into 
the  signal  by  means  of  a  second  lent  at  thown  in  Fig.  24;  at  pl«ae  B.  «  a 
suitable  tpatial  filtering  function  mi^  be  introduced  in  order  to  enhance  the 
signal- to-noite  ratio  for  a  desired  pattern.  For  example,  assume  that  a 
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transparency  T(x’,y') 

r'(x»  y‘)  f{x»,  y')  j 

ox  x(x,  y)  and  t  (x^ 


is  placed  at  plane 
the  image  produced 

y)  ;  i.e.. 


Pg  f  thus  producing  the  product 
at  plane  is  the  convolution 


f  (x,  y)  *  t  (x,  y)  - 


II f  (s,  n)  t  (x  -5 ,  y  -n)  d|  dn  (3k) 


It  is  therefore  possible  to  modify  the  spatial  frequency  content  of  a  light 
signal  distribution  in  a  manner  analogous  to  that  of  electrical  filters 
which  modify  the  tenq>oral  frequency  content  of  electrical  signals.  '  In 
particular,  using  this  technique,  one  can  utilize  optimum  smoothing  filters 
if  a  contour  obscured  by  other  detail  is  to  be  detected  in  its  original  form, 
or  matched  filters  if  recognition  of  the  pattern  is  sufficient;  in  the  latter 
case,  the  filter  acts  to  suppress  the  noise  and  to  concentrate  the  energy  of 
the  signal  into  a  small  area. 

To  illustrate  the  operation  of  a  matched  filter,  consider  again  the  case 
of  a  signal  1  (x,  y)  having  Fourier  transform  (x*,  y')  ;  and  assume  that 
the  signal  is  superimposed  on  a  noise  background  n  (x,  y)  .  The  matched 
filter  is  a  transparency  having  the  following  transmittance  function: 


T  (x»,  y')  -  K  E^*(x*,  y*)  (55) 

At  the  focal  plane  Pg  ,  the  signal  consists  of  a  wavefront  whose  contour  is 
characterized  by  a  phase  distribution  ^  (x*  y* )  .  By  application  of  the 
matched  filter  (55)  t  the  phase  is  compensated  for,  resulting  in  a  constant 
value;  i.e.,  in  a  plane  wavefront.  This  wave  is  then  focussed  by  lens  I2 
into  a  point  in  plane  ^  P^  ,  having  large  intensity  and  located  at  the  coordl* 
nates  corresponding  to  the  original  position  of  the  signal  in  plane  Pj^  .  The 
matched  filter  also  operates  on  the  amplitude  of  the  Fourier  transform  of  the 
signal,  attenuating  the  noise  where  the  signal  is  weak  and  thus  enhancing  the 
resultant  signal-to-noisa  ratio. 

In  order  to  take  advantage  of  the  signal-processing  properties  of  opti¬ 
cal  systems,  it  is  necessary  to  use  monochromatid,  coherent  light  illumina¬ 
tion.  The  realization  of  eonqplex  matched  filters  poses  a  difficulty  because 
such  filters  are  generally  obtained  with  photographic  film  transparencies  in 
which  the  transmittance  can  only  be  positive  or  zero,  and  the  phase  cannot  be 
controlled  simply*  (Attempts  to  control  the  phose  by  variation  of  the  thick¬ 
ness  have  proven  awkward. )  A  much  simpler  design  of  complex  filters  has  been 
developed  recent ly^^«  ^3  , 

Assume  that  a  transparency  with  the  complex  transmittance  function 

»  (**»  y*)  -  I  T  (x»,  y*)  I  ^ 
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must  be  built*  One  makes  a  mask  having  the  real  transmittance  function 


Tp  +  I  T  (x',  y*)  I  cos  (x^  x*  +  0  (s*,  y»)  )  (57) 

This  function  consists  of  a  constant  bias  term  To  (which  is  necessary 
because  the  transparency  must  remain  positive  everyv/here) ;  and  of  two  com> 
plex  terms^  respectively 


j 

I 

) 

i 


1 

i 

r 

] 


1 

2 


1 

2 


T  (x»,  y*) 


T  (x‘,  y*) 


exp  j  (x^  x*  +  0  (x*,  y*)  )  “ 

I  f  (x»,  y')  *^0 

exp  -J  (x^  X*  +  0  (x*,  y*)  )  - 

5  T*  (x*,  y*)  e  o 


(58) 


placing  the  abo\e  transparency  at  focal  plane  Pg 
tlon  of  the  signal  distribution  function  ?  (x,  y; 
of  (57)  ;  i.e., 


one  obtains  the  convolu- 
with  each  of  the  terms 


To  ^  (x#  y)  +  t  (x  -  \  i  y)*  1  (x  -  Xp,  y)  +  t#  (-x^  +  X,  -y)*. 

(59) 

•f  (-X  +  x^,  -y) 


In  order  to  illustrate  the  derivation  of  the  results  of  Equ,  59.» 
recall  that  multiplication  of  a  real  function  by  an  exponential  corresponds 
to  a  translation  of  variable  in  the  complex  frequency  domain;  i.e.,  if  !^(x'} 
Fou^r  transform  of  7(x)  ,  then  P(x'  -  a)  is  the  Fourier  transform 
)  ;  simllarlyi  ^(x*  +  a)  is  the  Fourier  transform  of  l*(x)* 


Thus,  the  derivation  of  Eq.  59  ^  convolution  of  the  signal  function 
with  each  term  of  the  transmittance  (57)  clear;  the  first  term  of  Eq.  59 
is  proportional  to  the  original  signal,  attenuated  by  the  constant  factor 
To  t  the  second  term  is  the  desired  matched  filter  output,  displaced  by  a 
distance  x^  froa  the  signal  of  the  first  ter%  and  the  third  term  is  still 
another  filter  output  obtained  with  a  filter  function  which  is  the  conjugate 
of  that  of  the  second  term,  displaced  a  distance  -x  .  The  displacement 
x^  is  purposely  introduced  in  order  to  avoid  overlapping  of  the  three  terms 
oi  Eq.  59  to  permit  their  separation, 

( 


is  the 

•• 
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Recapitulating,  we  have  shown  that  a  transparency  having  a  real  trans¬ 
mittance  function  type  (Eq.  57)  used  to  provide  a  complex  matched 

filter  operation,  l^e  preparation  of  the  transparency  (Eq,  57) 
by  various  methods;  an  elegant  approach  (Ref. 12)  is  one  based  on  t)  3  use  of 
a  modified  Mach-Zender  interferometer  (Fig.  25).  The  function  E(x,y)  , 
whose  Fourier  transform  is  Equ.  56,  is  placed  in  one  beam  of  the  :^nterfer- 
ometer,  so  that  the  corresnondlng  Fourier  transform  produced  by  lens  is 
displaced  at  its  back  focal  plane  outside  the  Interferometer,  Indicating 
with  r(x,y)  the  light  distribution  from  a  reference  beam,  whose  Fourier 
transform  is  R(x*,y*)  -  R  exp  j  (x^  x*  +  y  y*)  ,  one  has  the  folloi7ing 
resultant  light  distributiSn  at  the  output  p¥ane  : 


j0(x*,y») 

(60) 


At  plane  Pp,|  a  photographic  plate  with  suitable  exposure  slope  is  placed; 
this  gives  a jtranspareney  proportional  to  the  energy  output;  i.e,, 

2  '  2  2 

G  (x*,y*)  =  lR(x*,y*)  +  T(x*,y*)  |  *=|  R{x*y’)  j  +  |T(x',y*j  + 

+  2  I  R(x',y'')|‘|  T(x*,y‘)  |  cos  x*  +  y^  y*  +  0(x’,y')  ) 

(61) 


Substituting  the  expression  of  R(x',y*)  ,  one  finds* 

G(x’,y*)  ■  R^^  +  j  f(x‘,y’)  2R^|  T(x',y*)|  cos  (x^^  x*  + 

(62) 

+  y*  +  ^(x*jy*)  ) 


Thus,  the  required  transfer  function  type  (Eq.  57)  is  obtained. 


-21* 


(ONFIDENTIAl 


CONFIDENTIAL 


6.  DETECTION  OF  SIGfiATURE  BY  ELECTRONIC  FILTERING  PROCESSES  OF  PATTERN 

RECOGNITION 

When  viewing  a  picture,  the  human  eye  Is  found  to  be  more  sensitive  to 
the  detection  of  edges  (l*e.,  of  abrupt  rather  than  of  continuous  charges  of 
brightness).  This  property  has  been  utilized  In  the  study  of  television 
bandwidth  reduction  systems^^  .  With  reference  to  Fig*  26,  the  video  signal 
obtained  by  scanning  a  picture  Is  spilt  Into  two  paths— 'One  routed  through  a 
low-pass  filter  of  frequency  response  L(j(o),  and  the  other  one  routed 
through  a  differentiator.  The  former  portion  Is  limited  to  the  Information 
corresponding  to  the  slow  varying  signal  Intensity,  while  the  latter  provides 
a  large  output  in  correspondence  of  abrupt  signal  variations,  by  forming  the 
derivative  ^  ,  where  x  Is  the  scanning  variable  proportional  to  time. 
Thus,  the  latter  channel  contains  predominantly  edge  Information.  At  the 
receiving  end  of  the  TV  transmission  system,  the  edge- Information  signal  is 
passed  through  a  filter  whose  transfer  function  is; 


H  (J<d) 


1  ■  L 

^<D 


(63) 


The  output  of  the  latter  Is  added  to  the  slowly  varying  portion  of  the  infor¬ 
mation  and  reproduces  the  original  picture. 


A  generalization  of  the  method  of  edge  detection  has  been  proposed 
recently^^  with  application  to  two-dimension  processing.  The  method  Is 
illustrated  In  Fig.  27  and  utilizes  a  gradient  operator,  consisting  of  a  pair 
of  two-dimensional  filters  Hx  and  Bg  i  which  synthesize  the  high-frequency 
part  of  the  picture.  Using  appropriate  transfer  functions  H^(Ju,  Jv)  and 
Hp(jn,  Jv)  at  the  receiving  end,  and  recombining  the  outputs  as  In  the  case 
of  Fig.  26,  one  reconstitutes  the  original  picture. 


These  methods  have  the  advantage  of  requiring  a  smaller  bandwidth  than 
that  of  conventional  TV  transmissions;  thus,  bandwidth  reductions  as  high  as 
7:1  and  1^:1  have  been  demonstrated,  Exanq>le8  of  application  are  shown  In 
Pig.  28^®  . 


While  the  above-described  method  was  developed  predominantly  for  the  pui^ 
pose  of  bandwidth  reduction,  it  can  be  readily  modified  for  the  recognition 
of  V-shaped  patterns  such  as  are  of  interest  in  the  detection  of  optical  sig¬ 
natures  of  submarines. 
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7.  EXPERIMENTAL  WORK  IN  CONNECTION  WITH  THE  INVESTIGATION  OF  OPTICAL  AND 

ELECTRONIC  FILTER  TECHNIQUES  OF  PATTERN  RECOGNITION 

For  the  study  of  optical  matched  filter  and  electronic  edge-enhancement 
filter  techniques,  an  experimental  system  consisting  of  a  laser  Illuminator, 
beam  scanner,  receiver-collector  and  photodetector,  electronic  rate  counter, 
staircase  cathode-ray  tube,  display  signal  generator,  memory-type  oscillo¬ 
scope  display  has  been  built. 

The  block  dl  agram  of  the  laboratory  model  is  shown  In  Fig.  29*  It  con¬ 
sists  of  a  CW  visible  laser  (HeNe,  X  «>  6328^  ,  I*  ®  5  ®W)  whose  collimated 
output  is  scanned  mechanically  by  means  of  a  pair  or  reflecting  prisms  with 
axes  perpendicular  to  each  other  (Fig.  30);  prisms  are  driven  x^ith  step- 
type  motors  which  are  actuated  by  pulse  trains.  The  repetition  rate  of  each 
pulse  train  may  be  varied  at  present  up  to  values  of  45  pulses /sec  for  the 
faster  unit  (line  scanner),  and  the  repetition  rates  of  the  tv;o  pulse  trains 
are  ah^ays  maintained  with  a  constant  ratio  15:1  .  In  the  present  unit,  each 
pulse  produces  a  rotation  of  3*  of  the  axis  of  the  prism;  thus,  after  I5 
pulses,  a  total  rotation  of  45°  is  obtained,  and  a  nex^  face  of  the  prism  is 
presented  to  the  laser  beam.  While  one  prism  rotates  45°  >  the  other  one 
rotates  3°  j  so  on.  Thus,  one  prism  acts  as  a  line-stepper,  and  the 
other  one  as  a  frame-stepper.  Increasing  the  number  of  faces  of  the  prisms 
and  reducing  the  rotation  per  step,  any  angle  of  scan  in  the  vertical  and  in 
the  horizontal  direction  can  be  obtained. 

The  pulse  trains  are  generated  by  a  clock  consisting  respectively  of  an 
astable  silicon-controlled  rectifier  circuit,  oscillating  at  the  faster  rate, 
and  by  a  monostable  silicon-controlled  rectifier  circuit,  which  is  triggered 
periodically  by  one  of  the  pulses  of  the  first  train,  maintaining  a  fixed  ratio 
of  15:1  between  the  two  trains. 

The  clock  generator  is  also  utilized  to  produce  a  precision  stair-case 
voltage  for  the  synchronous ' control  of  the  two-dimensional  scanning  of  the 
readout  scope.  The  latter  circuit,  which  was  developed  along  principles  be¬ 
lieved  to  be  novel,  is  Illustrated  in  Fig.  31*  Two  stair-case  waveforms  are 
obtained  respectively  for  the  horizontal  and  for  the  vertical  deflection  con¬ 
trol.  Each  waveform  consists  of  I5  levels  or  steps,  and  is  obtained  by  recourse 
to  a  nonlinear  network  consisting  of  diode-resistance  circuits  in  which  the 
biases  of  a  number  of  diodes  are  progressively  changed  from  "off”  to  "on",  or 
vice  versa,  by  means  of  pulse-controlled  SCR  pulse  generators.  With  reference 
to  the  schematic  of  Fig.  &  dc  voltage  is  applied  across  the  divider  from 
A  to  B  and  from  A'  to  B*  .  The  divider  is  tapped  at  I5  points  equally  spaced 
in  resistance  values. 

Initially,  when  power  is  first  turned  on,  the  voltage  at  Points  1  through 
15  are  at  zero  potential,  because  the  SCR's  are  not  conducting.  This  is  true 
for  both  horizontal  and  vertical  circuits.  Diodes  Dl  through  DI5,  and  Dl* 
through  Dl5'i  sre  forward-biased,  and  therefore  Diodes  dl  through  dl5,  and  dl* 
through  dl5',  are  reverse-biased;  as  a  result,  the  output  is  0  volts  for  both 
horizontal  and  vertical  outputs. 
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The  voltages  at  the  divider  taps  will  appear  one  at  a  time  across  the 
load,  respectively  H  and  V  terminals,  when. the  corresponding  diode  of  the 
Dl-to-Dl5  group  or  Dl*-to-Dl5*  group  hccomes  reverse-blased.  Since  the 
voltages  at  the  control  points  of  these  diodes  are  applied  sequentially,  the 
output  becomes  a  stair-case  function. 

A  ring-counter  circuit  is  used  to  apply  pulses  sequentially  to  the  con¬ 
trol  terminals  1  through  1^.  The  design  of  this  circuit  Is  based  on  the  use 
of  cathode  coupling:  it  will  be  noted  that,  with  this  type  of  coupling,  the 
"on"  SCR  does  not"turn  off"  until  the  next  SCR  is  "on";  this  feature  is 
essential  in  order  to  prevent  occurrence  of  "drop"  in  the  output  stair-case 
voltage  waveform. 

Tlie  clock  switches  the  ring  counter  from  one  stage  to  another,  and  con¬ 
sequently  switches  the  output  to  the  next  step.  The  clock  pulses  have  no 
effect  on  the  ring  counter  until  both  "set"  switches  are  momentarily  closed. 

Photographs  of  the  stair-case  voltage  and  of  the  corresponding  two- 
dimensional  scope  display  are  shown  in  Figs.  ^2  and  33  respectively. 

As  seen  from  the  block  diagram  of  Fig.  2%  the  received  laser  signal  is 
first  applied  to  an  optical  signal-processing  unit  and,  subsequently,  is 
detected  with  a  photomultiplier  (surface  Si  or  S20).  The  output  of  the  latter 
is  proportional  to  the  reflected  light  signal  and  may  be  applied  directly 
(Switch  S  closed)  to  the  z-axis  of  the  display-scope,  thus  providing  an  opti¬ 
cal  representation  of  the  reflection  from  the  sea  surface.  An  electronic 
picture-processing  unit  is  expected  to  be  interposed  between  photomultiplier 
and  display  scope  to  provide  pattern  enhancement  by  edge  detection.  In  the 
present  set-up, the  optical  processing  unit  and  the  electronic  processing  unit 
have  not  yet  been  developed.  This  mode  of  operation  of  the  Laser  Detection 
system  results  in  a  television- type  picture  presentation  and  may  be  associated 
with  a  fast  scanning  of  the  sea-surface  area  under  observation.  Hake  patterns 
are  extracted  from  the  noisq  backgroiind  represented  by  the  wind-generated  waves 
and  are  viewed  directly. 

An  additional  mode  of  operation  of  the  Laser  Detection  system  is  shown 
in  Fig*  31  io  obtained  by  opening  Switch  S;  in  this  case,  a  Rate  Counter 
unit  is  interposed  between  photomultiplier  output  and  electronic  signal- 
processing  unit,  to  provide  a  signal  proportional  to  the  rate  of  oscillation 
of  the  waves  rather  than  to  the  size  of  their  reflecting  facet.  The  Rate 
Counter  has  already  been  built;  it  consists  of  a  pulse-shaping  circuit,  in 
which  the  received  pulses  are  converted  into  rectangular  pulses  of  uniform 
amplitude  and  duration,  and  of  an  integrator  circuit  whose  output  voltage  is 
proportional  to  the  rate  of  oscillation  of  the  waves.  In  this  mode  of  opera¬ 
tion,  a  pattern  different  from  the  'Srake"  is  obtained  on  the  scope,  namely,  a 
pattern  representing  loci  of  constant  oscillation  rate;  since  the  mean  fre¬ 
quency  of  sea  waves  is  of  the  order  of  10*20  per  sec  for  capillary  waves  .and 
of  the  order  of  1-2  per  sec  for  gravity  waves,  the scanning  speed  of  the  read¬ 
out  laser  beam  must  be  adapted  accordingly,  with  a  value  lower  than  fbr  wake 
pick-up. 
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Photographs  of  the  equipment  built  and  described  above  are  shown  in  Fig.  3^ 
and  Pig. 35;  these  Include  the  Laser^  the  Scanner  and  a  small  water  tank, 
and  the  Rack-mounted  Receiving  unit  and  oscilloscope. 

In  Its  present  form,  the  Optical  Signature  Detection  system  provides  a 
15  X  15  matrix-type  display  which  is  intensity-modulated  as  a  function  of 
the  laser  reflection  from  the  water  surface.  An  example  of  pattern  taken 
with  an  optical  collector  of  small  aperture  is  given  in  Pig.  36;  only  a  few 
illuminated  spots  appear  on  the  scope  because,  due  to  the  absence  of  diffu¬ 
sion  at  the  reflection  from  the  surface  of  the  water,  the  returned  beam 
remains  collimated  and  misses  the  optical  collector  at  a  number  of  points. 

In  order  to  apply  the  techniques  of  optical  matched  filtering,  a  suit¬ 
able  transparency  must  be  built;  as  previously  indicated,  the  aim  of  the 
filtering  process  is  the  optlmli^ation  of  the  peak  signal  energy  with  respect 
to  mean-square  noise  energy.  If  one  indicates  with  ?(x,y)  the  signal  and 
with  n(x,y)  the  noise,  and  assumes  that  the  latter  is  a  homogeneous  iso¬ 
tropic  random  process  with  spectral  density  N(x*,y')  ,  the  solution  of  the 
problem  is  obtained  by  recourse  to  a  filter  having  the  following  transmittance 
function^^  ; 


H(x*  y') 


(64) 


where  J*(x*y*)  is  the  complex  conjugate  of  the  signal  spectrum  and  k  is 
a  constant.  One  realizes  the  above  transfer  function  by  means  of  two  trans¬ 
parencies  placed  in  contact  with  each  other;  the  first  transparency  repre¬ 
sents  F*(x',y*)  and  is  made  by  recourse  to  the  Mach-Zender  interferometer, 
previously  illustrated  (Fig.  25).  The  second  transparency  represents 
l/N(x',y')  and  is  made  from  a  photograph  of  the  distribution  n(x, y)  ,  first 
by  producing  the  Fourier  transform  N(x*,y')  and  then  recording  it  on  a  suit¬ 
able  film.  If  D  is  the  density  of  the  film,  the  specular  a;.iplitude  trans¬ 
mission  is  I  ■  exp  (-0/2)  f  The  relationship  between  density  D  and  expo¬ 
sure  E  of  the  film  is  given  by  the  following  relationship: 


D  ■  y  (log  E  -  log  Ep  )  (65) 

where  y'  ood  are  characteristics  of  the  film  and  of  its  development 
process;  therefore,  om  has: 


I 


(66) 


Kaklng  E  proportional  to  N(x',y*)  and  ■  2  ,  the  specular  trans¬ 
mission  T  is  proportional  to  l^N(x',y')  •  ' 
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The  transparency  proportional  to  F*(x',y*)  may  be  made  taking  as 
reference  photographs  of  typical  vake  patterns  in  the  absence  of  noise. 

It  is  of  interest  to  recall  some  results  of  pattern  recognition  reported 
in  the  literature^®  ,  Several  geometric  patterns  are  shc«m  in  Fig,  37 J  those 
with  symmetry  (such  as  the  rectangles)  have  a  real  Fourier  transform,  while 
those  without  symmetry  (such  as  the  letter  I.)  have  a  complex  Fourier  trans-  . 
form.  The  complex  Fourier  transform  of  the  letter  L  is  shown  in  Fig.  38; 
this  was  taken  by  means  of  a  Mach-Zender  interferometer  and  consists  of  vari¬ 
ous  components--!, e,,  a  center  component  corresponding  to  the  first  tv70  terms 
of  Eq.  61,  and  two  side  components  corresponding  to  the  two  exponential  com¬ 
ponents  of  the  last  term  of  Eq,  6l  (as  shovm  in  connection  with  Eq's,  57 
58,)  Finally,  the  cross-correlation  of  the  L  mask  with  the  input  function  is 
shown  in  Fig,  39  appears  as  a  brilliant  symmetrical  signal^'  , 
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8.  CONCLUSIONS  .  , 

The  problem  of  the  detection  of  the  optical  signature  of  submarines 
may  be  approached  successfully  by  means  of  laser  techniques.  The  general 
method  discussed  In  this  report  is  based  on  the  use  of  an  unmodulated  laser 
beam  which  illuminates  the  sea  surface  from  the  air  and  is  scattered  with 
variable  Intensity,  optical  frequency,  angular  direction,  polarization, 
depending  upon  fii  e  state  of  the  sea  surface.  In  particular,  the  optical 
frequency  variations  depend  on  the  relative  velocity  of  the  sea  waves  and 
may  be  measured  by  maans  of  doppler  techniques;  unfortunately,  the  appli¬ 
cation  of  these  is  extremely  difficult  and  Impractical  at  optical  frequency, 
due  to  the  fact  that  the  heterodyne  detector  acts  as  an  antenna  whose  aper¬ 
ture  A  and  angle  of  arrival  of  received  rays  ft  must  satisfy  the  condi¬ 
tion  ^  i  where  X  Is  the  optical  wavelength, 

A  much  mere  practical  method  of  detection  of  the  optical  signature  is 
one  based  on  the  detection  of  the  modulation  of  the  angular  direction  of 
scatter  of  the  laser  beam,  which  varies  periodically  with  the  oscillations 
of  the  sea  waves,  and  on  Che  realization  of  a  TV-lllce  representation  of  the 
sea  surface,  with  the  utilization  of  optical  matched  filters  for  the  enhance¬ 
ment  of  the  slgnal-to-nolse  ratio. 
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PIG.  3  -  STATISTICAL  DISTRIBUTION  OF  WAVE 
SLOPES  (Ref.  1) 
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EFFECT  OF  BEAM  WIDTH  OH  DOPPLER  FRE 
qUEHCY  VALUE 
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PIG.  10  -  BLOCK  DIAGRAM  OP  OPTICAL  LASER  RADAR 
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TIG.  11  -  ANGULAK  DEVIATION  OF  EEFLECTEO  RAT 


FIG.  12  -  ANALYSIS  OF  ANGULAR  DEVIATION  OF  REFLECTED 
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PIG,  13  -  SCHEMATIC  OF  HETERODYNE  SYSTEM  WITH  FRE¬ 
QUENCY  TRANSUTION 
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FIG.  Ih  -  OF  PRACTICAL  RCl^.CDYKE  0^  HETERCLYIJE  SYSTEM 
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?2C.  Xfi  -  a)  PLOT  OF  AJiPLITUDB  0?  D8TBCTED  SIGNAL 
VERSUS  bI>^/2 

b)  TYPICAL  ANGULAH  SENSITIVITY  OP  HETERO- 
DY-KE  DETECTION 
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FIG.  17  -  EXPERDEKTAL  DOPPLER  lASER  RADAR  SYSTEM 
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FIG.  18  -  TYPICAL  FM  SPECTAUM  CORRESPONDING 
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KOmATTON  0?  LOCAL  OSCILLATOR 
AT  100  kc,  AtJD  MODITLATIOS  OF 
BACK  SIDEBAND  C(X'<?OK£NT  AT  10  kc 
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ri6,  eo  -  ENLARGED  VIEW  OP  THE  SPECTRUM 
SUEROOSDTHC  ONE  0?  THE  SIDEBAND 
(XD^PCfNENTS  OF  Fig,  1% 
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no.  fel  -  SPKCTRUM  RK3KLTING  FROM  MODUUTION 
OF  THE  LOCAL  OSCXLUTOR  AT  100  U 
AW)  MOOUUTIOW  OF  THE  TAROBT  AT  5  kc. 
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PIG.  S2  -  EXPERI>'ENTAL  SET-UP  FOR  STUDY  OF  REFLECTTO:« 
FROM  WATER  SURFACE 
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FIG.  23  -  FOURIER  TRANSFORMATION  OF  SIGNAL  LIGHT 

DISTRIBUTION 


C 


PIG.  2k  -  TWO-LENS  SYSTE’'  FOR  FILTERING  AND  REPRODUCTION 
OP  SIGNAL  LIGHT  DISTRIBUTION 
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FIG.  25  -  preparation  op  COMPLEX  KATCHEO  FILTER  TRANS¬ 
PARENCY 
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PIG.  26  -  EDGE- ENHANCEMENT  SYSTEM  (l-dlmenslon) 
(Reference  l6) 
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FIG.  -  TWO-DIMENSIONAL  EDGE- ENHANCEMENT  SYSTEM 
(Reference  16) 


CONFIDENTIAL 


MG.  28  -  EXA>!PtE  OF  TWO*Dn!ENSIOM;\L  EDGE 
EKHAJJCEMEirr  PROCESS  (Ref.  I3) 
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FI6.  31  •  STAIRCASE  VOLTAGE  GEITERATOR  AND  RATE  CLOCK 
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PIG*  3^  ”  LASEUj  MI5C[iAI<ICAt  SCAWJSR,  RK8CH  KAXER  TAiJK  »ITIi  WAVE 
K4KER.  Aliio,  A  Photosiultipll&r  And  &  Fresnel  T,eQe 
Collector  are  shown. 


PIC.  35 
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KACK-KOUNTED  BQUIPlSttr  (RATE  CtOCK,  STAXRCASB 
CSKBRATOEf  AMPLIFIER,  RASE  COUNTER,  POWER  SUPPLY) 
AND  OXSPUY  SCOPE. 
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riG.  36  -  EXAMPIBS  OP  PATTERN  OISPIAY 
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ric,  ^  •  CCRttLBX  rOORTER  TRAKSPCWf  OF  L- 
PATTERN  (Rtf.  11?) 
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